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Abstract 

A combinatorial interpretation is provided for the moments of characteristic 
polynomials of random unitary matrices. This leads to a rather unexpected conse¬ 
quence of the Keating and Snaith conjecture: the moments of | £(1/2 + it) | turn 
out to be connected with some increasing subsequences problems (such as the last 
passage percolation problem). 


1 Introduction 

Keating and Snaith JlJ] have proposed to model the limiting distribution of the Riemann 
zeros using the characteristic polynomials of unitary random matrices U 

Z(U,9) = det(I - e~ ie U) (1) 

In particular they deal with a conjecture for moments of |<£(1/2 + it)\ which states that 
the limit 

T 

'»=& (bih :? / | « i/2+a ) |2 ’” dt < 2 > 

exists and is equal to a product of two factors, f(m) and a(m), i. e. 

I m = a(m)f(m ) (3) 

The first factor a{m) is the zeta-function specific part, 


a(m) = [IP 

p 


y£M r(fc + m ) \ 2 

p to V k[ r M ) P 


( 4 ) 
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where the product is taken over prime numbers p. As for the second factor f(m), Keating 
and Snaith |l| have hypothesized that it is the random matrix (universal) part and may 
be represented as 

f(m) = lim N-™ 2 (\Z(U,9)\ 2m ) u{N) (5) 

JN—>00 

Here the average is over Circular Unitary Ensemble (CUE) of random unitary matrices 
NxN. 

In this letter I provide a combinatorial interpretation for the moments of characteristic 
polynomials, (| Z(U, 0)\ 2m )u(N)- I then relate these moments with two-rowed lexicographic 
arrays which are generalizations of permutations and words in combinatorics (for basic 
information about permutations, words and lexicographic arrays see, for example, a book 
by Fulton j2]). The combinatorial interpretation of moments of the characteristic poly¬ 
nomials enables an explicit formula to be obtained for the total number of two-rowed 
lexicographic arrays constructed from the letters of an alphabet of m symbols with the 
increasing subsequences of the length at most N. 

This combinatorial interpretation is, in fact, a natural consequence of a profound re¬ 
lation between Random Matrix Theory and Robinson-Schensted-Knuth type problems 
discovered recently by Gessel 0, Balk, Deift and Johansson 0 and Rains 0 (see also 
Aldous and Diaconis jSJ and the references therein ). In particular, it follows that cer¬ 
tain expectation values over CUE appear in the theory of last passage percolation. The 
moments of characteristic polynomials , (\Z(U, 0)\ 2m )u(N), may also be considered in this 
context. The reasoning along those lines enables the random matrix part f(m) defined 
by equations (@)-(@) to be related with the weakly-right/weakly-up lattice model of the 
last passage percolation. 

2 Increasing subsequences of lexicographic arrays and 

(\z{u,ef m )u ( N) 

Let us recall a generalization of the Cauchy identity (see Gessel ||, Balk and Rains fjj, 
Tracy and Widom || ): 

S A (£l, 6, • • • , Cm) SA (hi, V2, = 

N> Ai...>A m >0 

(detl n (.I + m p + ^t/t))) ( 6 ) 

' V>i =1 / ' U(N) 

Here s\ (£i, £ 2 , ■ ■ ■, Cm) denotes the Schur polynomial of m variables. Take £i = £2 — ■ ■ ■ — 
Cm = e~ lG and 771 = ?/ 2 = ■ ■ ■ = Vm = in equation (| 6 ). We will use the formula for the 
Schur polynomials of identical variables, 

s\(x, x ,..., x) = x Xl+X2+ - +Xm d x (m) (7) 
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where the coefficient d\(m ) is calculated by the formula 


d\(m) 


n (Xi-i-Xj+j) 

l<i<j<m 

0! 1! 2!... (m - 1)! 


( 8 ) 


From equations (|§|)-(j5[) we conclude that the moments of characteristic polynomials, 
(\Z(U, 0)\ 2m )u(N)i may be represented as sums over partitions. Importantly, these sums 
should be taken only over partitions of the length less than m and with the first row less 
than N. An explicit expression for the m LI1 moment of the characteristic polynomial is 


+oo 

(\z(u,o)\ 2m ) uim = Y, 


K =0 


Y d l ( m ) 

Ah K, Ai <iV 


(9) 


where A h K means that the set (Ai, A 2 ,..., A m ) is a partition of K, i.e. Ai + A 2 + ... + 
A m = K. Representation (|9]) enables a combinatorial interpretation to be provide for 

<1 z(u,e)\ 2m ) u{N) . 

Consider semi-standard Young tableaux constructed from K boxes with at most m 
and N boxes in hrst columns and first rows respectively. The total number of pairs 
of such tableaux is then given by the sum X) d|(m). (If is a well-known fact ( 

Ah K, Ai <N 

see Fulton ||2[) that the coefficient d\(m) defined by equation (g) gives the number of 
semi-standard Young tableaux of the shape A whose entries should be taken from the 
set [1,2 ,... ,m\). It is this sum which appears in expression (|9J) for the moments of the 
characteristic polynomials. In what follows we will apply the Robinson-Schensted-Knuth 
correspondence (see Fulton |2[]) between two-rowed lexicographic arrays and pairs of semi¬ 
standard Young tableaux. 

By definition, a two-rowed array of the size K is an object of the form 


A 


m _ 

2 ,K ~ 


u 1 U 2 ... U K 
V\ V 2 ... v K 


( 10 ) 


where the letters rq, u 2 , ■. ., Uk and v\, v 2 ,..., Vk belong to an alphabet (any ordered 
set) of m different letters. The following two conditions ensure that the array A™ K 
corresponds to a pair of semi-standard Young tableaux (each of K boxes with entries 
taken from 1 , 2 , ,m): 

U\ < u 2 < ... < uk ( 11 ) 

Vi < v 2 < ... < v K (12) 

The arrays A' 2 K for which conditions (|TT| ) and (|T^ ) are satished are called lexicographic 
arrays. 

Let us define a (weakly) increasing subsequence of the lexicographic array A™ K as 
follows: 


1. 


The element of the array A™ K with the number i\ is less than that with the number 
* 2 , (“[[) < if Yi < V i2 and u n < u i2 . 
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2. A subsequence of s elements of the array A™ K such that 



ii<i 2 < ... <i k 


will be called a weakly increasing subsequence of the length s. 


(13) 


Let l(A? K ) be the length of the maximal (weakly) increasing subsequence of the two- 
rowed lexicographic array A™ K . Denote by Rf nN the number of two-rowed lexicographic 
arrays of the size K constructed from the alphabet of m letters and including weakly 
increasing subsequences of the length at most N, i.e. 


R 


K 

m,N 


of A™ k with 1{A™ k ) < N 


(14) 


The number R m ^ defined by the sum over size K, 


+oo 

Rm,N = Rm,N (15) 

K =0 

gives the number of arrays of an arbitrary size that include weakly increasing subsequences 
of N elements or less. By the Robinson-Schensted-Knuth correspondence, the number of 
ordered pairs of semi-standard Young tableaux with m and N boxes in the first column 
and the first row respectively is equal to the number R m ,N of lexicographic arrays. 

This observation enables R m ,N to be expressed as an average over CUE 


+oo 

Rm,N ^ ^ 
K= 0 


H d l ( m ) 

Al -K, Ai <N 


(\Z(U,0)\ 2m ) utm 


(16) 


Here we have used the representation of the moments {\Z(U, 0)\ 2m )u(N) as the sums over 
partitions, see equation (|^). We thus conclude that the moments of characteristic poly¬ 
nomials, (\Z(U, 0)\ 2m )u(N), have a clear combinatorial interpretation. They are equal to 
the number of two-rowed lexicographic arrays constructed from an alphabet of m symbols 
with the weakly increasing subsequences of the N elements at most. 

ft is now possible to obtain an explicit formula for the number R m ,N of lexicographic 
arrays. Once we apply the result of Keating and Snaith[l] : 


N 

(\ziu,e)\ 2m ) um = W 


3 =1 


r (j) r (j + 2m) 

[r(j +m)f 


(17) 


equality (0) between the number of arrays R m ,N and the moments of characteristic 
polynomials gives 


N 

Rm,N 

i=i 


r (j) r (i + 2y ») 

[r(i + m)] 2 


(18) 


In order to illustrate how formula (|18|) works consider the following example. Take 
an alphabet consisting of two letters a and b, i.e. = (a, b), m = 2. It may give 
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rise to both lexicographic and non-lexicographic arrays. A typical lexicographic array 
compounded from the letters of the alphabet (a, b ) is 

f a a b b \ 

\a b b b ) 

We observe that conditions (|TT|), ( |T2|) are satisfied by this particular array. An example 
of a two-rowed non-lexicographic array is given below: 

( a a b b \ 

\a b a b ) 

In this array the second letter of the word abab, b , is larger than the third letter of this 
word a. Thus, condition ([12|) is not satisfied resulting in a non-lexicographic array. 

Let us list all possible two-rowed lexicographic arrays from the alphabet (a, b) with 
the weakly increasing subsequences of the length at most N = 2: 



There are 20 arrays with the required properties. It can be easily verified that formula 
(|TS|) obtained above gives precisely this number, i.e. 

Rm,N = ft r |rr+ + W ) = 20 (™ = 2,N = 2) (19) 

j =1 [r(j+m)] 2 

3 Moments of characteristic polynomials in the last 
passage percolation theory 

Let us now turn to an interesting interpretation of lexicographic arrays in the framework 
of the last passage percolation theory. (See, for example, Baik j9j] where different last 
passage percolation models are described and the relations with different expectation 
values over CUE are explained). It is a consequence of the relation between the moments 
of characteristic polynomials and lexicographic arrays that (\Z(U, 0)\ 2m )u(N) appear also 
in the last passage percolation problems. 

Consider the lexicographic array A™ K of the form ( 0 ) where the letters Ui,Vi]i = 
1,2 ,,K take values in the set of integers = 1,2 ,,m. This lexicographic array 
corresponds to a matrix with non-negative integer entries. The (i,j) entry of this matrix 
is the number of times the vector (J.j occurs in the array A™ K . For example, the array 

( 1 1 1 2 2 2 3 3 3 3 \ 
[1221111233] 
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corresponds to the matrix 


/ 1 2 0 \ 

3 0 0 

VI 1 2 ) 

Let e [1, 2, m] be a planar array of non- negative integers. Consider weakly- 

up/weakly-right paths 7 r' s , (■ ik,jk)k=v suc h that R — R < ■ ■ ■ < k and ji < J 2 < • • • < ji 
(This model has discussed at length in Baik |Dj][|rOJ). Let (1,1) /* ( m,m ) denote the set 
of up/right paths from (1,1) to (m, m) in the planar array X(i,j). The entries of the 
matrix X(i,j), Xij , may be considered as the passage times at knots (i,j). Each planar 
array of non-negative integers X(i,j ) can be assigned the last passage percolation time 
T(m,m ) for travel from (1,1) to ( m,m ). T(m,m) is defined explicitly by the following 
expression: 

T(m, m) — max V (20) 

So far we have used the Robinson-Schensted-Knuth correspondence between pairs of semi¬ 
standard Young diagrams and lexicographic arrays. Turning to the last passage perco¬ 
lation problems Robinson-Schensted-Knuth correspondence is a correspondence between 
planar arrays X(i,j ) of non-negative integers and pairs of semi-standard Young tableaux. 
In this case the last passage percolation time T(m,m ) should be considered instead of 
the length of the longest increasing subsequence. The number of arrays R m ,N defined in 
the previous section by equations (B and (pT|) is replaced by the number of all possi¬ 
ble X(i,j) of size mxm with the last passage percolation time T(m,m ) less than N. 
Applying the results of the previous section we find that 

It of X(i,j) with T(m,m) < N = (\Z(U, 0)\- in ) u( ^ N) (21) 

With the help of the Keating and Snaith formula (equation (|I7j)) an explicit expression is 
obtained for the number of planar arrays X(i,j) of the size m x m with the last passage 
percolation time bounded by N: 


ft of X ( i , j ) with T(m, m) 


N 


<^=n 

1=1 


£ (i)r(i + 2 ra) 

[r (j + m)} 2 


( 22 ) 


4 Combinatorial consequences of the Keating and 
Snaith conjecture 

The above results lead to an interesting interpretation of the Keating and Snaith con¬ 
jecture. Once m is a non-negative integer, the random matrix (universal) factor f(m) 
relates the moments I m of zeta-function with increasing subsequences. On assumption 
that the Keating and Snaith conjecture is valid, that is equation (||) holds, it follows from 
the results of the previous sections that 

OO 

E It Of A™ k with 1{A£ k ) < N 

K =0 _ 

jy m 2 



Im , ■ 

——- = Inn 

a[m) N—>oo 
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In other words, the ratio defined by equations (|2|)-(f|) is asymptotically equal to the 
number of lexicographic arrays constructed from an alphabet of m letters with increasing 
subsequences of at most N elements divided by N m . 

As a particular case consider the weakly-up/weakly-right last passage percolation 
model. Here the ratio is equal to the whole number of m x m planar arrays X(i,j) 

of non-negative integers with last percolation time at most N divided by N m : 

Im = lim tt of X(i,j) with T(m, m) < N 
a(m ) N™ 2 ^ ' 

Apart from the combinatorial interpretation of the Keating and Snaith conjecture, the 
above considerations raise a question of how the moments of the Riemann zeta function 
are related to the representation theory of finite groups. 
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